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We study the long-time evolution of the ion temperature in an expanding ultracold neutral plasma using 
spatially resolved, laser-induced-fluorescence spectroscopy. Adiabatic cooling reduces the ion temperature by 
an order of magnitude during the plasma expansion, to temperatures as low as 0.2 K. Cooling is limited by 
heat exchange between ions and the much hotter electrons. We also present evidence for an additional heating 
mechanism and discuss possible sources. Data are described by a model of the plasma evolution, including the 
effects of ion-electron heat exchange. We show that for appropriate initial conditions, the degree of Coulomb 
coupling of ions in the plasma increases during expansion. 


I. INTRODUCTION 


Ultracold neutral plasmas (UNPs), formed by pho- 
toionizing laser-cooled atoms near threshold, have ion 
temperatures around 1K and tunable electron temper¬ 
atures from IK to 1000Ki One of the most important 
topics studied in UNPs is the physics of strongly coupled 
plasmas, which have a ratio of nearest neighbor potential 
energy to kinetic energy larger than one. This is quanti¬ 
fied by the Coulomb coupling parameter 
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where a ws = (47rn/3) is the average interparticle 
spacing (i.e. Wigner-Seitz radius) and T is the temper¬ 
ature of the particle species of interest. Ions in UNPs 
have T,; « 2 — 4, and they have been used to study the 
equilibration of strongly coupled Coulomb systems after a 
potential quench^ 8 . and collision rates beyond the regime 
of validity of Landau-Spitzer theory^ Strong coupling 
is of interest in many different plasma environments^ 
and the equilibration dynamics observed in UNPs in par¬ 
ticular has been connected to laser-produced plasmas 
formed from clusters and solids^i and proposed as a pos¬ 
sible limitation in the brightness of photoemitted electron 
beams— 

A central focus of research on UNPs is the creation 
of more strongly coupled plasmas to study collisions, 
transport, collective modes, and many other phenomena 
across a wider variety of coupling strengths into the liq¬ 
uid regime. Proposals to increase T^ by modifying the 
plasma after its creation are sensitive to the long-term 
evolution of the ion temperature. Here, we report mea¬ 
surement and modelling of the ion-temperature evolution 
of expanding UNPs, which has never been experimentally 
studied before because of the difficulty in distinguishing 
the small thermal motions of ions from the large, hydro- 
dynamic ion expansion velocity— - — Our results show 
significant ion adiabatic cooling, with ion temperatures 
decreasing up to an order of magnitude. We also observe 
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significant contributions to the ion temperature from ion- 
electron thermalization, as well as an additional source of 
ion heating we ascribe to deviations of the plasma density 
from an ideal spherical Gaussian. 

This study adds to prior studies of fundamental plasma 
properties conducted using UNPs, such as the plasma 
creation process,2il i 13 ’ 14 collective modes of ion s 16 : 17 and 
electrons , 13 i 18 i 19 and formation and ionization of Ryd¬ 
berg atoms in the plasma.— - — 

There is long-standing interest in the problem of a 
plasma expanding into vacuum, which typically dom¬ 
inates the dynamics of plasmas created with pulsed 
lasers,— such as in experiments pursuing inertial confine¬ 
ment fusion^ x-ray lasers^ or the production of ener¬ 
getic (>MeV) ions through irradiation of solids; 26 ’ 27 thin 
foilsand clusters— 

For ultracold neutral plasmas, an early hydrodynamic 
model and numerical study, which included the ef¬ 
fects of three-body recombination and other inelastic 
processes^ accurately described experimental observa¬ 
tions of the ion density evolution— An adiabatic, self¬ 
similar solution of the Vlasov equation s 15 ’ 31- — describes 
the expansion dynamics well across a wide range of initial 
conditions, and this is the starting point for understand¬ 
ing the experiments described here. This model predicts 
adiabatic cooling of electrons, which has been confirmed 
through the expansion dynamics ; 15 ’ 35 electron loss^li and 
the TBR rate— The model predicts adiabatic cooling of 
the ions as well, and it was pointed out that this should 
result in an increase of correlations to Tj « 10 before 
the density and collision rate drop so much that correla¬ 
tions freeze out ] 38 ’ 39 Small effects of ion correlations on 
the expansion have been discussed in conjunction with a 
comprehensive model^ 8 Our work provides the first test 
of these predictions of the evolution of the ion tempera¬ 
ture. 

Adiabatic cooling has been observed in trapped non¬ 
neutral plasmas when reducing confinement for both pure 
ion^H and antiproton plasmas— “2, It is important in the 
interstellar medium and solar winder— and it has been 
discussed for ions in a plasma created by laser irradiation 
of a solid— >22 
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II. EXPERIMENTAL METHODS 

UNPs are created by photoionizing laser-cooled 88 Sr 
atoms in a magneto-optical trap.— This results in a spher¬ 
ical Gaussian density distribution (initial e -1 / 2 density 
radius = op ~ 1 mm to 2 mm) of a few hundred mil¬ 
lion ions (and electrons) with a peak density of up to 
5 x 10 16 m -3 . 

Spatially and temporally resolved laser induced fluo¬ 
rescence spectroscopy on the ion’s principal transition at 
422 nm permits in situ probing of the local kinetic energy 
and density of the ions / 1 5,48 We excite fluorescence in a 
1 mm thick sheet that bisects the plasma. Fluorescence 
emitted perpendicular to the sheet is collected by a 1:1 
optical relay and 4x objective/ocular magnification stage 
that images onto an intensified charge coupled device 
(ICCD), with a pixel size of 13 pm. This allows regional 
analysis of small volumes of the plasma with roughly con¬ 
stant density and bulk expansion velocity— The total 
system resolution is mainly limited by the ICCD, which 
suffers from blooming during the amplification process. 
To minimize this effect, we experimentally obtained the 
point spread function of the ICCD^ and deconvolve each 
image in post analysis — - 

Creation and probing of the plasma is repeated with 
a 10 Hz repetition rate, and the LIF-laser frequency is 
scanned to build up an excitation spectrum. The spec¬ 
trum of the region of interest is fit to a Voigt profile. 
We assume a Lorentzian-component width dominated by 
the laser and natural linewidths (6 MHz, 20 MHz, respec¬ 
tively) and the Doppler-broadened Gaussian width is a 
fit parameter related to a local ion temperature.— The 
signal is proportional to density and can be calibrated 
using absorption imaging^ with a resulting uncertainty 
in density of approximately ±20%. 

To minimize the contamination of ion temperature 
measurements by expansion velocity^ temperatures are 
measured for regions that are only one pixel wide along 
the LIF beam axis and 1 mm perpendicular to it. To im¬ 
prove statistics, we analyze 79 regions covering a 1 mm 2 
area, and compute an average temperature from the mea¬ 
surements. Error bars on temperature measurements are 
statistical uncertainty in the resulting mean. The averag¬ 
ing over multiple regions introduces a maximum relative 
density variation of e -1 / 8 = 0.88 from the center region 
to the outermost region, for op =lmm. All densities we 
quote are an average of the regional densities, > 95 % of 
the peak density. 


III. OVERVIEW OF ULTRACOLD PLASMA 
EVOLUTION TIME SCALES 

The plasma is unconfined and undergoes a complex 
evolution as it expands into the surrounding vacuumji 
with dynamics that can be divided into three distinct 
time scales as shown in Figure [I] 

Global thermal equilibrium for electrons is established 
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FIG. 1. Three main stages of UNP evolution. Electrons 
reach global thermal equilibrium (GTE) on a timescale ap- 
proximat ely equal to the inverse electron plasma frequency, 
ujpe = y/eom £ /n e e 2 . Ions come to local thermal equilibrium 
(LTE) on a timescale of the inverse ion plasma frequency. 
Hydrodynamic effects such as plasma expansion occur on the 
longest timescale. 


on a very fast timescale compared to the physics of in¬ 
terest in this paper at a temperature determined byi^S 

3ksT e 

—= hun- E t . (2) 

Here, tvjji is the combined energy of the photoionizing 
photons and is the ionization energy of a ground state 
atom. We can create UNPs with initial electron temper¬ 
atures of 1 K to 1000 K, however we typically stay above 
40 K to avoid three body recombination effects, which 
become important when T e > 0.1.— We ionize a max¬ 
imum of about 25% of the atoms for T e < 150 K, but 
this approaches 100% near an auto-ionizing resonance 
(Te ~ 430 K) — All but a few percent of the electrons 
are trapped by the space-charge field of the ions— 

Next, the ions come to local thermal equilibrium on 
a time scale comparable to the inverse of their fu n- 
damental oscillation frequency, w” 1 = y/eomi/riie 2 ~ 
1 ps. Initial ion velocities and positions are inherited 
from the laser-cooled atoms, resulting in low kinetic en¬ 
ergy and high potential energy because the spatial dis¬ 
tribution is completely random. Coulomb energy is con¬ 
verted into ion kinetic energy in a process called disorder- 
induced heating (DIH) or correlation heating,—^ which 
yields a density-dependent, equilibrium ion temperature 
(To/ij oc n -1 / 3 ~ IK). Disorder-induced heating limits 
the ion’s Coulomb coupling parameter to 2 < Fj < 4, 
with the variation determined by electron screening— 
Strictly speaking ions never achieve global thermal equi¬ 
librium in our experiments because of the long timescale 
for thermal transport— 

In the next stage of evolution, the plasma expands 
into the surrounding vacuum in a hydrodynamic fashion, 
driven by the electron thermal pressure. For a quasi¬ 
neutral, perfect Gaussian density distribution, assuming 
global thermal equilibrium of ions and electrons and neg¬ 
ligible inelastic collision processes and electron-ion ther- 
malization, the expansion is described by a self-similar 
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expansion 

a 2 (t)=a 2 0 (l+t 2 /r 2 exp ), (3) 

T < ,e(t) = T <1 e(0)/(l + t a /O J (4) 

where the expansion timescale, r exp = 
^/mitTg/A:B[T e (0) + T((0)] « 10 ps, is set by initial 
size and temperatures after electrons and ions reach 
global and local equilibrium respectively ! 15 i 32 The 
electron temperature dominates, which allows one to 
indirectly measure T e ( 0) by fitting the evolution of 
the cloud size. Equations [3] and [4] also show that 
T e ,icr 2 = constant , reflecting the adiabatic cooling of 
both species as the cloud expands. We will discuss 
below how these equations must be modified to account 
for heat exchange between electrons and ions, but they 
provide good intuition for plasma dynamics. 




time (/is) time/(27r/oj p j) 


FIG. 2. Early time measurements of ion kinetic energy widths 
for various density UNPs with T e = 430 K (a) Disorder in¬ 
duced heating temperatures vary with density as shown, (b) 
Scaling temperature by the predicted DIH magnitude and 
time by the ion plasma oscillation period results in close to 
universal behavior. The solid lines are results from molecular 
dynamics simulation of an equilibrating plasma which are fit 
to experimental data to determine ion density. 


IV. INITIAL ION EQUILIBRATION: 
DISORDER-INDUCED HEATING, KINETIC-ENERGY 
OSCILLATIONS, AND DENSITY CALIBRATION 

The initial equilibration of the ions has been exten¬ 
sively discussed previously^—— but we describe it here 
in order to give a complete picture of the temperature 
evolution and because early-time heating dynamics pro¬ 
vide the most accurate determination of the plasma den¬ 
sity. 

Figure [5] shows the increase of ion kinetic energy due to 
initial disorder (DIH) and the subsequent kinetic energy 
oscillations (KEOs) at close to twice the ion plasma fre¬ 
quency ujpi that result from the evolution of spatial corre¬ 
lations in the equilibrating strongly coupled plasma. The 
equilibrium temperature from DIH for a homogeneous 


plasma, is set by density (oc n 1 / 3 ) and electron tempera¬ 
ture (via electron-screening). It can be calculated fromS 


Tdih = 


2 e 2 
3 kjs 47r£o aws 
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using tabulated molecular dynamics simulation results^ 
for U, the excess particle energy in units of e 2 /Atieq »«»■ 
k = a ws /\D is the screening parameter for electron De¬ 
bye length Ad = (eofcsTe/nee 2 ) 1 / 2 . A closed-form de¬ 
scription of KEO dynamics does not exist, and it can 
only be predicted with molecular dynamics simulations. 
The oscillation frequency has been shown by particle- 
in-cell Yukawa simulations^ tree-code algorithms^ and 
full molecular dynamics simulations^ to agree well with 
u>pi, although electron screening softens the ion-ion in¬ 
teraction and slows the oscillations—^ This deviation is 
small but observable for our conditions. As seen in Fig. 
[2Kb), when time is scaled by 27r/w p i, and ion temperature 
is scaled by Tdih, the various plasma exhibit similar be¬ 
havior, and oscillations damp to a scaled temperature 
close to unity. 

All presented temperatures are fit parameters assum¬ 
ing a Voigt spectral profile with Maxwellian velocity dis¬ 
tribution. It has been shownfZ however, that the velocity 
distribution shows a non-Maxwellian, high-velocity tail 
that relaxes on the same timescale as the oscillations, so 
Voigt fits will underestimate the ion RMS kinetic energy 
until the ions have achieved local thermal equilibrium. 
Therefore, to interpret our early-time measurements we 
compare them to a library of ion velocity distributions 
from molecular dynamics simulations of a homogeneous 
plasma— 2 We use the numerical data to produce simu¬ 
lated LIF spectra that are fit with the same procedure 
as experimental data (Fig. [2] (b)). Assuming an electron 
temperature given by the photoionizing-laser photon en¬ 
ergy we can allow density to be a fit parameter, and the 
fit is highly constrained because the KEO frequency and 
the temperature both depend on the density. Good confi¬ 
dence in the measured density and Tdih is important for 
accurate simulation of the full ion temperature evolution 
as discussed below. Note that the measured temperature 
only approaches Tdih as the oscillations damp. 

Fits of experimental data to numerical results are pre¬ 
sented as insets throughout the paper. A full study of 
this density-fitting method will be given elsewhere— - It 
is more accurate than using LIF intensity calibrated with 
absorption imaging^ but the two methods agree well, 
yielding an absolute uncertainty in density on order of 
± 10 %. 


V. ION ADIABATIC COOLING 

On a hydrodynamic timescale of ~ 10 ps the plasma 
cloud will expand due to the thermal pressure of the elec¬ 
trons, leading to adiabatic cooling of electrons and ions. 
When ion temperatures are measured for longer evolu- 
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tion times we clearly see cooling of the ions by up to an 
order of magnitude, as demonstrated in Figs. [3] and SI 
To establish that we are observing adiabatic cooling of 
the ions, we check the scaling predicted by Eq. [4] Ac¬ 
cording to this simple model, the only two parameters 
determining the ion temperature are the initial temper¬ 
ature, set by DIH, and the expansion time T exp . Figure 
[3] (b) shows the evolution of the ion temperature scaled 
by those two parameters for various initial electron tem¬ 
peratures holding plasma size and initial peak density 
fixed. Figure [4] (b) shows a similar study varying the 
plasma size for similar initial electron temperature and 
ion density. In both cases, the scaled data approximately 
collapse onto universal curves, showing that adiabatic 
cooling, described by Eq. [4] dominates the ion temper¬ 
ature evolution. We observe temperatures as low as 0.2 
K, which are the lowest temperatures ever measured for 
ions in local thermal equilibrium in an ultracold neutral 
plasma. 


VI. MODELLING UNP DYNAMICS 

To more accurately model the data, we use a hybrid ap¬ 
proach first developed in Ref— i that combines a hydrody¬ 
namic treatment of the evolution of the plasma size and 
electron and ion temperatures with terms derived from a 
kinetic and numerical treatment of the effects of correla¬ 
tions. The model takes the form of differential equations 
describing the evolution of the Gaussian size parameter 
u, the electron and ion temperatures, and the expansion 
parameter 7 . The hydrodynamic expansion velocity for 
the ions at position r with respect to plasma center, is 
v(r, t) = 7(t)r. The differential form of the equations 
allows inclusion of new terms to treat electron-ion ther- 
malization. 


( 6 ) 

dl (t) k B T e (t) + Uii(t)/3 ^ 

—5T“ = - T77\ -7 (t) {<) 

at rrii(T z {t) 

® j 1 = - 27 ((W(t) -§ (~tmm++ 

2^ lei (t)T e (t) 

nrii 

( 8 ) 

= -2 7 (i)T e (t) - 2^7 ei(t)T e (t) (9) 

at rrii 


where Unit) is the average ion-ion correlation energy per 
particle in the plasma, which can be related to the spatial 
pair correlation function of the ionsji Uu is negative in 
UNPs. 

The electron-ion equilibration terms in Eqs. [ 8 ] and [51 
which contain 7 e j will be discussed below. When they are 
omitted and the correlation energy is neglected (Uu = 0 ), 
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FIG. 3. Ion temperature evolution for UNPs with similar ini¬ 
tial plasma size (oo = 1 mm) and comparable ion densities 
(rii(0) = 2 x 10 15 m~ 3 to 4 x 10 15 m -3 ) and different initial 
electron temperatures given in the legend, (a) Lower ini¬ 
tial electron temperature produces stronger electron shield¬ 
ing and lower temperature after DIH. (b) Scaling the time 
and temperature axes yields universal curves that only vary 
with electron screening. The comparison of experimental data 
with molecular dynamics simulations (solid lines) is used to 
determine the local ion density, (c) Lower electron tempera¬ 
ture yields slower expansion and ion cooling, (d) Scaling the 
temperature by Tdih and time by the expansion time T exp 
approximately collapses the data onto a universal curve, con¬ 
firming that adiabatic cooling dominates the long-time evolu¬ 
tion. The line is the prediction of the ion temperature model 
('Eas.Hfl- nUl) including correlation effects and cloud expansion 
only,— omitting electron-ion thermalization 


Eqs. [G][9] can be solved exactly by the analytic solution 
presented in Eqs. [3] and [4] 

The inclusion of the correlation energy in Eqs. [G] - [9] 
allows an approximate treatment of correlation effects 
such as disorder-induced heating. Uu can only be calcu¬ 
lated with a full molecular dynamics description of the 
plasmai It’s evolution can be approximated, however, 
as2£ 


dUu(t) _ Uu(t) - U u , eq (t) 

dt ~ T corr (i) ’ 1 j 

where the timescale for relaxation of the correlation en¬ 
ergy to its equilibrium value, Uu, eq r^ is taken as the 
inverse-ion plasma frequency T corr = y/mi£o/e 2 rii. 
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FIG. 4. Long time ion temperature evolution for UNPs with 
similar initial electron temperature (T e (0) = 430 K) and com¬ 
parable ion densities yet different initial plasma sizes (details 
in the legend), (a) Plasma size does not effect local DIH dy¬ 
namics. (b) Comparison with molecular dynamics simulations 
is used to determine local ion density, (c) Larger initial cloud 
size yields slower expansion and ion cooling, (d) Scaling the 
temperature by Tdih and time by the expansion time r exp 
approximately collapses the data onto a universal curve. 


At plasma formation, Ua is set to zero, but its magni¬ 
tude increases rapidly, raising the ion temperature close 
to the equilibrium temperature Tdih on a timescale of 
T CO rr- The ion temperature after establishment of local 
equilibrium is not put into the model as an ad hoc pa¬ 
rameter; it emerges naturally from the dynamics. This 
treatment also includes the retardation of the expan¬ 
sion due to decreasing magnitude of Ua with decreasing 
density, which is a small effect for our conditions. The 
ion KEOs observed experimentally during the disorder- 
induced-heating phase are not described at this level of 
treatment. The slow approach to Tdih during the damp¬ 
ing of the KEOs mentioned in Sec. IIVI is also not de¬ 
scribed. 

Inherent in the traditional interpretation of Eq.[8l how¬ 
ever, is an assumption of global thermal equilibrium for 
the ions that is not well-justified. The equilibrium corre¬ 
lation energy depends on the local density, and disorder- 
induced heating only leads to local thermal equilibrium. 
The characteristic timescale for global ion thermal equi¬ 
libration Ttherm is set by the thermal conductivity for 
the strongly coupled ions, which can be approximated 


as A ~ nkBOJpia^,,,^- Temperature equilibration between 
the ions in two regions of volume F, interface area A, and 
separation d is ~ cyVd/XA , where cy = inks/ 1 } is the 
ion heat capacity per volume. Taking all length scales as 
the plasma size <r, this yields T t herm ~ (J 2 /wpja^ s , which 
is on the order of 10~ 3 s for a typical UNP - much longer 
than other dynamic timescales for the plasma. 

This creates a significant complication for a rigorous 
description of the plasma, but we take advantage of the 
extremely long timescale for ion global equilibration to 
make a significant approximation to describe our data. 
We restrict our analysis to a central region of the plasma 
in which the average density is 96% of the peak density. 
We consider Eq. [S] as describing an effective tempera¬ 
ture for local thermal equilibribrium within this region. 
While previous treatments^ have approximated Uu teq (t) 
as an average over the equilibrium correlation energy in 
the entire plasma, we take it as the correlation energy 
for the measured plasma density, which is the average 
density of our region of interest. Several additional fac¬ 
tors make Eq. [8] a reasonable approximation for the local 
ion conditions. The actual error introduced with this ap¬ 
proximation is small because the temperature after local 
equilibration ~ Tdih varies slowly with density. Also, 
the self-similar plasma expansion (Eqs. [3]and[]]) leads to 
the same relative change in volume and adiabatic cooling 
in all regions of the plasma. Finally, the plasma expan¬ 
sion, which controls the adiabatic cooling dynamics, is 
dominated by the electron temperature for our condi¬ 
tions and is relatively insensitive to the ion temperature. 

With this interpretation, the equilibrium value of the 
correlation energy ( Uu, eq < 0) is estimated as^ 

Uu,eq{Ti,n) = k B TiT 3/2 {Ti,n) 

( A\ _i_ ^ (H) 

\ + r(Tj, nj 1 + r(Tj, n) j ' 

where X) and n are the instantaneous local ion temper¬ 
ature and density, and A\ = —0.9052, A 2 = 0.6322, and 
A 3 = — V^3/2 — A\/ \J A‘i. 

Eas. iGl fTUl neglecting electron-ion equilibration qualita¬ 
tively predict the ion temperature evolution, as shown 
by rescaling data by Tdih and r exp and overlaying the 
results of the model (solid lines in Figs. 0 (b) and0](b)). 
But quantitatively, the theory clearly diverges from the 
data after the DIH phase and underestimates the ion 
temperature during the expansion. This implies that the 
model misses a significant amount of energy that is trans¬ 
ferred to the ions during the evolution, and we now turn 
to a discussion of the sources of this energy. 


1. Electron-Ion Thermalization 

Collisional energy transfer between electrons and ions 
is typically neglected in theory and experiments on 
UNPsji Due to the large mass difference, complete ther¬ 
malization between the two would require 10 to 100 times 
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the expansion time ( T exp ) of an UNP. Thus, the energy 
lost by electrons due to collisions with ions is insignif¬ 
icant compared to the electron temperature,— and the 
expansion of the plasma and the electron temperature 
evolution can be described by Eqs. [3] and [4] 

However, if one is concerned with the ion temperature, 
electron-ion thermalization must be considered since the 
ions are typically orders of magnitude colder than elec¬ 
trons. This is particularly important for proposed ef¬ 
forts to increase Coulomb coupling by laser coolin g 63 i 64 
or circumventing DIH!! because a small transfer of elec¬ 
tron energy to thermal ion energy heats the ions signif¬ 
icantly. In general, thermal relaxation in a plasma has 
been a long-standing fundamental interest^— as well as 
a critical aspect of many areas of study such as inertial 
confinement fusion;' 3-7 ' warm dense matter^ and space 
plasmas.— 

The measurements presented here are the first to follow 
the ion temperature on the timescale of the expansion 
and observe the effects of electron-ion thermalization in 
an UNP. To describe the data, and help us identify the 
contribution of this heating mechanism, we start with the 
classic electron-ion collision frequency in a singly ionized 
plasmas, assuming rii = n e , nii » m e and T e >> Ti^ 

4\/27m^e 4 In [A] 

3(47t £ 0 ) 2 m 1 e /2 (k B T e ) 3 / 2 

= yAr 3 /2 Wpe l n[A ], (12) 

and the ion heating and electron cooling rates due to 
electron-ion collisions (EICs), 


dTj _ _dTe _ TOe 
dt dt rrii ^ el e 


(13) 


The argument of the Coulomb logarithm is A = 
1 / . The Coulomb coupling parameter for the 

electrons is never greater than 1 for our experiments, 
so strong coupling effects are smallas has been 
shown by comparison of Eq. [T5] with quantum T-matrix 
calculations!! in this regime of coupling. 

The heating and cooling terms are included in the evo¬ 
lution of ion and electron temperatures (Eqs. [8] and [9]). 
We use the measured, average plasma density for the re¬ 
gion of interest in these expressions. Following a similar 
line of reasoning as discussed in the context of the corre¬ 
lation energy, this provides a reasonable approximation 
of the evolution of the local ion temperature in the cen¬ 
tral region of the plasma because of the very slow global 
equilibration of the ions. 

With the inclusion of electron-ion thermalization in the 
model, it is possible to determine the significance of the 
effect for ion temperature evolution. Figure [5] shows data 
and simulation for various regimes of UNP parameters. 
The simulation results show that the heating from EICs is 
more significant for higher density, smaller electron tem¬ 
perature, and larger initial plasma size. These conditions 
increase the collision rate and the time (~ T exp ) before 


the density drops and collisions become negligible. In 
Fig. E5b), electron-ion thermalization contributes 0.5K 
and doubles the ion temperature at later times. For the 
opposite extreme of the accessible regime of these param¬ 
eters in UNPs, the effect is very small. 



FIG. 5. Comparison of data with models including and not in¬ 
cluding electron ion collisions. The lines are simulation results 
of the ion temperature model (Eqs. 161 101) including correla¬ 
tion effects and cloud expansion!! with (solid) and without 
(dotted) electron-ion thermalization. The insets show early 
time KEOs with time normalized by 2The solid lines 
in the insets are results from MD simulations fit to the data 
in order to determine experimental density. 


2. Additional Heating/Ion Acoustic Waves 

It is evident from Fig. [5] that the simulation still un¬ 
derestimates the measured ion temperature during the 
expansion by as much as several hundred millikelvin, and 
that the excess energy appears on the timescale of r exp . 
We observe that this effect is very sensitive to experimen¬ 
tal parameters such as laser alignment and beam quality. 
But in general the additional energy is more significant 
for larger density and electron temperature. 
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One possibility is that hydrodynamic expansion ve¬ 
locity is being misinterpreted as thermal velocity. The 
plasma is expanding with a velocity that varies in space, 
and LIF signal captured from a region with finite width 
Sx will reflect a spread of expansion velocities that can 
broaden the measured spectrum. For the self-similar ex¬ 
pansion (Eqs.[3]|3]), the velocity spread Sv = 7 (t)Sx peaks 
at t = r exp and then decreases. Assuming an ideal single 
pixel region of width Sx — 13 pm, Sv peaks at a value of 
lms -1 corresponding to a Doppler broadening of 2 MHz 
or 8 mK in temperature units. The actual value would in¬ 
crease if the expansion deviates from the ideal self-similar 
expansion or if the point spread function used for image 
deconvolution is inaccurate or not constant throughout 
the fluorescence volume. However, it is unlikely this ac¬ 
counts for a large fraction of the observed heating. 

Another possibility is that Eqs. [TS] and [13] underesti¬ 
mate the amount of ion heating due to electron-ion col¬ 
lisions. We do not attempt to use our data to extract 
the electron-ion thermalization rate because of the diffi¬ 
culty in separating this heat source from other possible 
mechanisms. While there are many assumptions that go 
into the derivation of Eqs. [T5] and [T3] it is unlikely that 
this can be the source of the anomalous heating because 
it appears to scale very differently than heating due to 
electron-ion collisions (Fig. [5|. 

We hypothesize that the largest contribution to the 
observed excess ion energy reflects ion acoustic waves 
(IAWs) that are excited during the ionization process 
due to deviations from a smooth Gaussian in the initial 
density distribution of the plasma. Previously we have 
used transmission masks on the ionizing beam to inten¬ 
tionally modulate the initial density distribution lead¬ 
ing to the excitation of IAWsi 16 ’ 85 ’ 86 However, similar 
excitations would be caused by any imperfection in the 
ionizing laser profiles and/or atom density distribution. 
Long-wavelength IAW excitations with a broad spread 
of wavevectors would increase the RMS velocity width 
of the ions. Short wavelength excitations should de¬ 
cay quickly through Landau damping, increasing the ion 
thermal energy. Both would increase the temperature as 
measured by our LIF probe. 

To assess the feasibility of this explanation, we esti¬ 
mate the energy per ion from a single IAW mode of fre¬ 
quency w, wavevector fc, and fractional density modula¬ 
tion — as^£ 

n 0 

A E iaw « imi (j — \ . (14) 

4 \k n 0 J 

In the long wavelength limit ui/k can be replaced with 
the sound speed y/ksTe/nii. If we assume all the IAW 
energy is transferred into thermal energy, we find that 

A T„^\T e (—\ . (15) 

6 \n 0 J 

For a ^ of only 5% and T e = 430 K this is nearly 200 mK 
of energy. In the short wavelength limit to/k is replaced 


by ojpi/k leading to a decreasing energy with wave vector, 


AT, 


1 ne 2 
6 £q ksk 2 



(16) 


This effect seems capable of providing the observed 
heating, and residues of the density-distribution fit to 
a 2D Gaussian are often on the order of 5 % on a ~ 
1 mm length scale. Similar depth ion holes have been 
shown to perturb the ion temperatures by several hun¬ 
dred millikelvim 17 ! 87 

Despite great effort, deviation of the density distri¬ 
bution from an ideal Gaussian remains a limitation for 
achieving lower ion temperatures. For all but the lowest 
atom densities, the magneto-optical trap for the atoms 
is optically thick to the laser-cooling light, this results 
in non-Gaussian atom distributions that are inherited by 
the plasma. Also, the output of the pulsed-dye laser that 
ionizes the atoms has spatial intensity modulations that 
are imprinted onto the plasma. Complete saturation of 
the ionizing process could minimize this source of heat 
but it would require larger pulse energies than available 
with our current system. 


VII. INCREASE IN T z WITH EXPANSION 

For the best chance of seeing increasing with 
expansion/^ minimal heating and fast expansion is re¬ 
quired. Fast expansion is achieved by small cloud size 
and high T e . Higher density doesn’t change the ini¬ 
tial coupling parameter, but it increases the DIH tem¬ 
perature - minimizing the relative decrease in com¬ 
ing from additional heating sources. Higher density also 
enhances EIC heating, but for sufficiently high electron 
temperatures, e.g. near the auto-ionizing resonance, the 
EIC cross section is so small that the highest achievable 
densities should not experience significant EIC heating. 
Higher density also allows data collection at later times. 
Therefore, at least for our accessible experimental param¬ 
eters, high T e , small cloud and high density are optimal 
for observing increasing Tj. 

Measurements at later times are limited by plummet¬ 
ing densities and signal-to-noise ratios, but for the best 
conditions, we see increase above 5, as shown in Fig. [G] 
EICs and other sources of heating prevent Tj from reach¬ 
ing as high as expected from Eqs. ICTTU1 neglecting these 
effects. The measured value of is calculated from ob¬ 
served ion temperature and density, and we are not able 
to determine if spatial correlations exist corresponding to 
equilibrium at this value. —- 


VIII. CONCLUSION 

We have used spatially resolved fluorescence spec¬ 
troscopy to measure ion temperature evolution during 
the expansion of an UNP for various initial plasma sizes, 
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FIG. 6. Increase in the ion Coulomb coupling parameter with 
plasma expansion for T e (0) = 430 K, no = 4 x 10 15 m“ 3 , and 
no = 1 mm. The line is the prediction of Eqs. I6IT01 


electron temperatures, and densities. We have observed 
several effects in UNPs for the first time, such as ion adi¬ 
abatic cooling and ion heating due to collisions with elec¬ 
trons. We have presented a numerical model describing 
the evolution of plasma size and electron and ion temper¬ 
atures including the effects of electron-ion collisions. We 
also presented evidence for additional ion heating, which 
we hypothesize is dominated by IAWs excited during 
plasma creation. We see increasing coupling for ions with 
time, but it is limited by these heating effects. Attempts 
to access stronger Coulomb coupling, such as through 
proposed schemes to laser cool ions in UNP o 63,64 ' 88- — or 
simply through plasma expansion as demonstrated here, 
should benefit from improved control of the initial den¬ 
sity distribution of the plasma. 
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